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Abstract 



We briefly review several models of heavy quarkonium production in hadronic col- 
lisions, and discuss the status of QCD factorization for these production models. 
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1 Introduction 



RHIC has produced sufficient evidence that a new state of hot and dense 
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matter of quarks and gluons, or the quark-gluon plasma (QGP), was formed 
in ultrarelativistic heavy ion collisions [1]. Anomalous suppression of heavy 
quarkonium, J/ip, was suggested as a good probe for discovering the QGP [2]. 
Because of the heavy mass, production of heavy quark pairs takes place at a 
very short time and is unlikely to interfere with the formation of QGP. On 
the other hand, the probability for the heavy quark pair to become a bound 
meson in the QGP could be suppressed due to the Debye screening from color 
charges in the medium. 

However, the suppression of heavy quarkonia in a dense medium also depends 
on how the bound state was formed. A quicker formation would give less time 
for the medium's color charges to break the coherence of the heavy quark pair, 
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and would require a higher color density and/or temperature to produce the 
suppression. It is therefore very important to ask how well we can calibrate the 
suppression in heavy ion collisions to extract useful information on medium 
properties. In this talk, I first briefly review several models for heavy quarko- 
nium production, and then, discuss the status of QCD factorization for these 
production models, and finally, give a brief summary and outlook. 



2 Production models 

Heavy quarkonium production has been the subject of a vast theoretical lit- 
erature and of intensive experimental study [3]. In order to produce a heavy 
quarkonium in hadronic collisions, the energy exchange in the collisions has 
to be larger than the invariant mass of the produced quark pair (> Iuiq with 
heavy quark mass mg). Since the binding energy of a heavy quarkonium of 
mass M is much less than heavy quark mass, (M 2 — Aitlq) / Aitlq <C 1, the 
transition from the pair to a meson is sensitive to soft physics. The quan- 
tum interference between the production of the heavy quark pairs and the 
transition process is powerly suppressed by the heavy quark mass, and the 
production rate for a heavy quarkonium state, H, up to corrections in powers 
of 1/rriQ, can be factorized as, 

o-a+b^h+x ~ / dT QQ a A+B^QQ[n]+x(T Q Q,m Q ) F Q Q [n] ^ H (Y q q) (1) 
n J 

with a sum over possible QQ[n] states and an integration over available QQ 
phase space oITqq. The nonperturbative transition probability F for a pair 
of off-shell heavy quark (tp) and antiquark (%) to a quarkonium state H is 
proportional to the Fourier transform of following matrix elements 

£<0|x f (yi) Kn HV2)\H + N)(H + N\^(y 2 ) K! n X {y,) |0> , (2) 

N 

where yiiyji) are coordinates, and IC n and JC' n are local combinations of color 
and spin matrices for the QQ state n. A proper insertion of Wilson lines 
to make the operators in Eq. (2) gauge invariant is implicit [4]. In Eq. (2), 
a A+B^QQ[n]+X) represents the production of a pair of on-shell heavy quarks 
and is calculable in perturbative QCD [5]. The debate on the production 
mechanism has being focused on the transition from the pair to the meson. 

Color-singlet model The color-singlet model (CSM) assumes that only a color 
singlet heavy quark pair with the right quantum number can become a quarko- 
nium of the same quantum number and the transition from the pair to a meson 
is given by the quarkonium wave function [6]. o"a+_b^qq+x' s dependence on 
relative momentum of the pair is neglected, and / dYqQ Fqq^ h in Eq. (1) is 
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Fig. 1. Charmonium cross section as a function of pt along with the CDF data 
points [8] and the theory curves from NRQCD model from Ref. [7] (a), and CEM 
from Ref. [10] (b). 

equal to the matrix element in Eq. (2) evaluated at yiiyi) = 0, which is pro- 
portional to the square of coordinate-space quarkonium wave function at the 
origin, \R H (0)\ 2 [3], and therefore, 



&A+B-+H+X ^ 'A+B-*QQ[H]+x( m Q) \Rh(0)\ 2 ■ 



(3) 



The same wave function appears in both production and decay, and the model 
provides absolutely normalized predictions. It works well for J/ip production in 
deep inelastic scattering, photon production, and some low energy experiments 
[3], but fails to predict the CDF data, see the dotted lines in Fig. 1(a) [7]. 

Color evaporation model The color evaporation model (CEM) assumes that 
all QQ pairs with invariant mass less than the threshold of producing a pair 
of open-flavor heavy mesons, regardless their color, spin, and invariant mass, 
have the same probability to become a quarkonium [9]. That is, the -Fqq 
in Eq. (1) is a constant for a given quarkonium state, H, and therefore, 

i-2Mq 

o-aIb^h+x ~ h / dm Q Q cr A+B ^ QQ+x (m Q Q) (4) 

J2m,Q 

with open- flavor heavy meson mass Mq and a constant /# [9] . With a proper 
choice of niq and the model gives a reasonable description of almost all 
data including the CDF data, as seen in Fig. 1(b). 

Nonrelativistic QCD model The Nonrelativistic QCD (NRQCD) model is 
based on the fact that the typical heavy quark rest-frame kinetic energy and 
binding energy, rriQV 2 , in a heavy quarkonium is much smaller than the heavy 
quark mass. The model separates the physics at scales of order mq and higher 
from the dynamics of the binding by using NRQCD, an effective field theory 
[11]. It provides a systematic prescription to calculate the physics at rrig order- 
by-order in powers of a s , and expands the nonperturbative dynamics in terms 
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Fig. 2. From Ref. [7], NRQCD predictions of charmonium polarizations are com- 
pared with the CDF data [12] . 

of local matrix elements in power series of heavy quark velocity v [7,11], 

NRQCD v"^ s _ //nH\ /r\ 

a A+B->H+X — 2-^ a A+B^QQ[ri\+X \ u n I > \°) 

n 

where the are NRQCD operators for the state H [11], 

Qf(0) = x t (0)X^^(0) (aW) ^ f (0)/C; X (0), (6) 

where o) H is the creation operator for H, x (VO are two component Dirac 
spinors, and K, n and /C^ are defined in Eq. (2) and can also involve covariant 
derivatives. At higher orders in v, the operator can have additional de- 
pendence on field strength as well as more fermion fields. The factorization 
in Eq. (5) could be understood from Eq. (1) by expanding the aA+B^QQ[n]+x 
at heavy quark relative momentum, q = (jpQ —pq)/2 = 0. The moments, 
J dTqQ q N Fqq^ h , lead to local matrix elements with high power of v. 

The NRQCD model allows every QQ[n] state to become a bound quarkonium, 
while the probability is determined by corresponding nonperturbative matrix 
elements (O^). Its octet contribution is the most important one for high pt 
quarkonium production at collider energies [7]. The NRQCD model has been 
most successful in interpreting data [3,7], as seen in Fig. 1(a). 

Quarkonium polarization and other models The key difference between the 
NRQCD model and the CEM is the prediction on quarkonium polarization. 
Once the matrix elements (O^) are determined, the NRQCD model can sys- 
tematically calculate the polarization of produced heavy quarkonia. We cannot 
calculate the polarization in CEM because of the uncontrolled radiation from 
the quark pair. The polarization of quarkonia at large px was considered a 
definite test of the NRQCD model [7]. But, the model has failed the test, as 
seen in Fig. 2, if the CDF data hold [12]. Several improved and new models 
have been proposed to address the issues of polarization [13]. 
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(a) (b) 

Fig. 3. Sketch for process A + B — > H(p)+X: (a) sample scattering amplitude with 
possible factorization breaking soft interactions indicated by dashed lines, and (b) 
factorization with Wilson lines indicated by thin double lines. 

3 Factorization 



Heavy quarkonium production in hadronic collisions involves both perturba- 
tive and nonperturbative scales. Nonperturbative physics appears not only in 
the transition from the heavy quark pair to a bound state but also in incom- 
ing hadron wave functions. A typical scattering amplitude for quarkonium 
production, as sketched in Fig. 3(a), can have soft and nonperturbative inter- 
actions between incoming hadrons as well as between the spectators and the 
formation process. These soft interactions may introduce process dependence 
to the nonperturbative matrix elements, and consequently, spoil the predictive 
powers of Eqs. (4) and (5). 

A proof of the factorization needs to: 1) show that the square of the scattering 
amplitude in Fig. 3(a), after summing over all amplitudes with the same ini- 
tial and final states, can be expressed as a convolution of the probabilities, as 
sketched in Fig. 3(b); each probability represents a square of sub-amplitudes 
and is evaluated at its own momentum scale(s); 2) show that the piece eval- 
uated at perturbative scale (s) is infrared safe and those evaluated at nonper- 
turbative scales are universal. 

As argued in Ref. [5], cross sections for producing on-shell heavy quark pairs 
can be computed in terms of QCD factorization. Therefore, the right-hand-side 
of the CEM formalism in Eq. (4) is perturbatively calculable for quarkonium 
production at a large pr, and the /# should be universal within the model be- 
cause soft interactions in Fig. 3(a) are suppressed by powers of l/pr- However, 
the factorization statement does not provide justification that the transition 
from a heavy quark pair to a quarkonium state is independent of the pair's 
invariant mass itiqq, spin, and other quantum numbers. When *C mg, the 
universality of /# may not be valid because of the spectator interactions. 

Fully convincing arguments have not yet been given for NRQCD factorization 
formalism in Eq. (5) [7,11]. Since the spectator interaction between the beam 



5 



jet and the jet of heavy quark pair should be suppressed by powers of 1/pt, one 
might expect the NRQCD factorization formalism to work at large pt- When 
Pt 3> ttiq, heavy quarkonium production is similar to the single light hadron 
production, and is dominated by parton fragmentation. The cross section is 
proportional to the universal parton-to-hadron fragmentation functions [4], 

<t a +b^h+x(pt) = ^A+B^i+x{PT/z^)®D H / i (z,m Q) p) + 0{m 2 H /p 2 T ) . (7) 

i 

Here, <g> represents a convolution in the momentum fraction z. The cross sec- 
tion aA+B^i+x includes all information on the incoming state, including con- 
volutions with parton distributions of hadrons A and B at factorization scale 
/i, as sketched in Fig. 3(b). As a necessary condition for NRQCD factorization 
in Eq. (5), the following factorization relation, 

D H/i (z,m Q ,n) = di^ Q Q [n] {z,fi,m Q ) (Off) , (8) 

n 

is required to be valid to all orders in a s and all powers in t>-expansion for 
all parton-to-quarkonium fragmentation functions [4]. In Eq. (8), g^qq[„] de- 
scribes the evolution of an off-shell parton into a heavy quark pair in state [n] , 
including logarithms of /i/mg, and should be infrared safe [4]. 

The factorization relation in Eq. (8) was tested up to next-to-next-to-leading 
order (NNLO) in a s at v 2 order in Ref. [4], as well as at finite v in Ref. [14]. 
Consider representative NNLO contributions to the fragmentation process of 
transforming a color octet heavy quark pair to a singlet, as sketched in Fig. 4. 
The individual classes of diagrams in Fig. 4(1) and (II), for which two gluons 
are exchanged between the quarks and the Wilson line, satisfy the infrared 
cancellation conjecture of Ref. [11], by summing over the possible cuts and 
connections to quark and antiquark lines, as do diagrams that have three 
gluon-eikonal vertices on the quark pair and one on the Wilson line [4]. For 
Fig. 4(111) type of diagrams, however, with a three-gluon interaction, this 
cancellation fails. By summing over all contributions to second order in the 
relative momentum q, it was found that the fragmentation function has a 
noncanceling real pole [4] 

ir i \v) = -a 2 s ±- £ v 2 . (9) 

From Fig. 4(111), this infrared divergence is not topologically factorizable 
and can not be absorbed into the matrix elements (Off). As demonstrated 
in Ref. [4], as a necessary condition for restoring the NRQCD factorization 
at NNLO, the conventional 0(v 2 ) octet NRQCD production matrix elements 
(Off) must be modified by incorporating Wilson lines that make them man- 
ifestly gauge invariant, so that the infrared divergence in Eq. (9) can be ab- 
sorbed into the gauge-completed matrix elements. 
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Fig. 4. Representative NNLO contributions to g — > QQ fragmentation in eikonal 
approximation, see Ref. [4] for the details. 

For quarkonium production, all heavy quark pairs with invariant mass less 
than a pair of open-flavor heavy mesons could become a bound quarkonium. 
Therefore, the effective velocity of heavy quark pair in the production could 
be much larger than that in decay. For charmonium production, charm quark 



velocity, v c ~ |<? c |A™c < y (4M|, — Am 2 c ) / (Aral) ~ 0.88, is not small, and 

therefore, the velocity expansion for charmonium production may not be a 
good approximation, unless one can identify and resum large contributions 
to all order in v or have a factorized formalism at finite v. It was found in 
Ref. [14] that with the gauge-completed matrix elements, infrared singularities 
in the fragmentation function for a color octet pair to a singlet at NNLO are 
consistent with NRQCD factorization to all orders in v or to a finite v, 



where f(v) = 2v/(l + v 2 ). The result in finite v is remarkably compact and 
intriguing, and should encourage further work on the factorization theorem. 



4 Summary and outlook 

Heavy quarkonium has two intrinsic scales and could be a good probe for the 
properties of the QGP or other dense medium. However, after more than 30 
years since the discovery of J/^, we still have not been able to fully understand 
the production mechanism of heavy quarkonium in high energy collisions, in 
particular, the transition from the heavy quark pair to a bound quarkonium. 
Although there are good reasons for each production model, none of the fac- 
torized production formalism, including that of the NRQCD model, has been 
proved theoretically. Further work on the factorization theorem is critical. 

The transition from the pair to a bound meson is sensitive to the QCD con- 
finement, and it is the dynamics of the transition that interacts and probes the 
properties of dense medium in heavy ion collisions. Nuclear matter could be an 
effective filter to distinguish the production mechanism [15]. Detailed study of 
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nuclear dependence of heavy quarkonium production in hadron-nucleus colli- 
sions should provide invaluable information on the formation of heavy quarko- 
nia in hadronic collisions. How good a probe the heavy quarkonium production 
is completely depends on how well we understand and are able to calibrate 
the production. 
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